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A systematic procedure to derive shell models for MHD turbulence is proposed. It takes into 
account the conservation of ideal quadratic invariants such as the total energy, the cross-helicity 
and the magnetic helicity as well as the conservation of the magnetic energy by the advection term 
in the induction equation. This approach also leads to simple expressions for the energy exchanges as 
well as to unambiguous definitions for the energy fiuxes. When applied to the existing shell models 
with nonlinear interactions limited to the nearest neighbour shells, this procedure reproduces well 
known models but suggests a reinterpretation of the energy fiuxes. 

PACS numbers: 



I. INTRODUCTION 

Understanding the existence and the dynamics of the 
magnetic field of the Earth, of the Sun and, in general, of 
other celestial bodies remains one of the most challenging 
problems of classical physics. Astronomical and geophys- 
ical observations have provided many insights into these 
phenomena [l|, 0, S]- Laboratory experiments 0,0] have 
confirmed that generation of magnetic field (dynamo) can 
take place under various circumstances, and lead to a 
variety of complex behaviours. However, analytical ap- 
proaches of this problem are extremely complicated while 
numerical efforts are limited to a range of parameter 
space that is often quite distant from the realistic sys- 
tems. For instance, in certain astrophysical bodies as 
well as in laboratory experiments, the kinematic viscos- 
ity V of the fiuid is six orders of magnitude smaller than 
its resistivity 77. The two dissipation processes therefore 
takes place at very different time scales. This property 
makes direct numerical simulation of dynamo intractable. 
Due to this reason we resort to simplified models. 

Shell models specifically belong to this class of simpli- 
fied approaches [0|. They have been constructed to de- 
scribe interactions among various scales without any ref- 
erence to the geometric structure of the problem. They 
were first introduced for fiuid _turbulence with the quite 

and have been extended 
|. In shell models, dras- 



successful GOY shell model Lz 
to MHD turbulence (9|,[l3,[ll 
tically reduced degrees of freedom (usually only one com- 
plex number) are used to describe the entire informa- 
tion provided by a shell of Fourier modes in wavenumber 
space. This approach reduces the description of turbu- 
lence from a partial differential equation to a reduced set 
of ordinary differential equations and provides a simpli- 
fied tool for studying the energy and helicity exchanges 
between different scales at a significantly reduced numer- 
ical cost. 

The present work aims at deriving the expressions for 
the energy fiuxes and the energy exchanges for magne- 
tohydrodynamic (MHD) turbulence in a systematic and 
consistent manner. Then we apply this scheme to study 
energy transfers in a shell model of MHD. This approach 



follows quite closely the previous efforts in which fluxes 
and energy exchanges have been identifled for the com- 
plete MHD equation. However, in these works [13, llJ|, 
energy exchanges between two degrees of freedom have 
been determined from the triadic interactions up to an 
indeterminate circulating energy transfer. The strategy 
adopted in the present paper is somewhat different. Here, 
we derive the energy transfer formulas from the energy 
equations by identifying the terms that participate in 
these transfers. This process also involves various sym- 
metries and conservation laws of the ideal (dissipation- 
less) equations. For example, the energy transfer from 
a magnetic fleld shell to another magnetic fleld shell is 
dictated by the convective term of the induction equa- 
tion that conserves the total magnetic energy. One of the 
main advantages of the present formalism is that we need 
not worry about the indeterminate circulating transfer 
appearing in the related past work by Dar et al. [13, llJ|- 



The dynamo process involves growth of magnetic en- 
ergy that is suppHed from the kinetic energy by the non- 
linear interactions. As we will show in the paper, a 
clear and unambiguous identiflcation of the various en- 
ergy fluxes and energy exchanges between the velocity 
and the magnetic fields is very important in the study of 
dynamo effects. This is one of the main motivations for 
the development of the present approach. The approach 
is also explicitly appHed in SectionJTVl to the derivation 
of the GOY shell model to MHD flir 



An outline of the paper is as follows: A general for- 
maHsm for expressing the various constraints satisfied by 
the nonlinearities in the shell models is discussed in Sec- 
tion [III It is shown in Section Hill that this formalism can 
be adapted nicely to the derivation of explicit expressions 
for the energy fluxes as well as for the shell-to-shell en- 
ergy exchanges in shell model. In Section IIV[ we apply 
the formalism to the GOY shell model for MHD turbu- 
lence [l3|, and study the energy energy fluxes for MHD 
turbulence. In Section V, we present our conclusions. 



II. SHELL MODELS OF MHD TURBULENCE 

Shell models were first introduced for fiuid turbulence 
(see for example 0, @, lla|)- They can be seen as a dras- 
tic simplification of the Navier-Stokes or the MHD equa- 
tions which, assuming periodic boundary conditions, are 
expressed in Fourier space as follows: 

-T— = nk(u, u) -nk(b,b) -z/fc^Uk + fk, (1) 

^ = nk(u,b)-nk(b,u)-77A:2bk. (2) 

where Uk and bk are the velocity and magnetic field 
Fourier modes respectively with wave vector k. The norm 
of this wave vector is A: = |k|. The viscosity v and the 
magnetic diffusivity ry are responsible for the dissipative 
effects in these equations while energy is injected through 
the forcing term fk. The nonlinear term is defined by 



nk(x,y)=zP(k) 



p+q=-k 



(3) 



where x and y can be either the velocity or the magnetic 
field. The tensor P is defined as 



P,-(k) = ^ 



fZj K j 



(4) 



It projects any field to its divergence-free part and it is 
used since only incompressible fiows are considered in this 
study (V-u = 0). In the velocity equation, the projection 
of the nonhnear terms using the tensor ^ replaces the 
introduction of the pressure term. In the magnetic field 
equation, the nonlinear terms are usually not projected to 
their divergence-free parts. Indeed, the non divergence- 
free parts of the two nonhnear terms cancel each other 
and the constraint V • b = is automatically satisfied. 
The writing of the nonlinear term in the magnetic field 
equation using the form ^ has been used to stress and 
explore the inner symmetries in the MHD equations. 

The incompressible MHD equations are known to con- 
serve the total energy, the cross helicity and the magnetic 
helicity. The conservation of these quantities plays a cen- 
tral role for the derivation of shell models. Similarly, the 
conservation of both the kinetic helicity and the kinetic 
energy in absence of magnetic field are used to simplify 
further the shell model for MHD. There is however an- 
other property that has not been exploited so far: The 
conservation of magnetic energy by the first nonlinear 
term in the magnetic field equation. Indeed, assuming 
periodic boundary conditions, it is easy to prove that 



Y^ nk(u,b) -bk^O. 



(5) 



The identification of a similar term in shell models for 
MHD will prove to be very useful in determining the en- 
ergy exchanges and the energy fiuxes in the shell model. 
The equations of the evolution of the variables in a 
shell model are designed to mimic as much as possible 



the MHD equations (fT]|2l) . In order to build the shell 
model using a systematic procedure, we first introduce 
the partition of the Fourier space into shells Si defined as 
the regions |k| S [h-i, ki] where ki = fco A*. In this def- 
inition, ko corresponds to the smallest wave vector. The 
number of shells is denoted by iV, so that the wave vec- 
tors larger than ko A^~^ are not included in the model. 
Any observable that would be represented in the origi- 
nal MHD equation by its Fourier modes Xk is described 
in the framework of the shell model by a vector of com- 
plex numbers noted X. Each component Xi of this vector 
summarises the information from all the modes Xk corre- 
sponding to the shell Si . It is also very useful to introduce 
the vector X; for which all components but the z-th are 
zero: 



X = (xi,X2, ■■■,Xn) e C^ , 
X, = (0,0,...,0,x„0,...,0)eC^ 



N 



X = 5]X,;, 



(6) 

(7) 

(8) 



where the expansion |[8l) is a direct consequence of the 
definition of X^ . 

In the following, the scalar product of two real fields 
will be needed for defining various quantities like kinetic 
and magnetic energies, cross hehcity, kinetic and mag- 
netic helicities. Using the Parseval's identity, the shell 
model version of this physical space scalar product is ex- 
pressed as follows: 



(X|Y) 



^ 1 

E^( 



XiVl+ViX*). 



(9) 



Due to the nonhnear evolution of the velocity and the 
magnetic field in the MHD equations, any attempt to 
design a mathematical procedure that would reduce the 
description of these fields to two vectors of complex num- 
bers U and B must lead to closure issues. In the deriva- 
tion of a shell model, the shell variables are usually not 
seen as projected versions of the original MHD variables 
and their evolution is not derived directly from the MHD 
equations ^l^ . The evolution equations for U and B are 
rather postulated a priori, but a number of constraints 
are imposed on the shell model. In this section, the 
models are build by imposing on the evolution equations 
for these vectors as many constraints as possible derived 
from conservation properties of each of the terms appear- 
ing in the original MHD equations. 

Property 1 ; the nonlinear term in the evolution equa- 
tion for \] is a sum of two quadratic terms; The first one 
depends on U only and conserves the kinetic energy £^ 
and the kinetic helicity Ti'' independently of the value of 
the field B; The second term depends on B only. 

Property 2 ; the nonlinear term in the evolution equa- 
tion for 'B is a sum of two bi-linear terms; The first one 
must conserve the magnetic energy 8^ independently of 
the value of the field U; 



Property 3 ; the full nonlinear expression in both the 
equations for U and B change sign under the exchange 
U^B; 

The dynamical system for the shell vectors can therefore 
be written: 

rfiU = Q(U,U)-Q(B,B)-iyD(U) + F, (10) 
dtB = W(U,B)- W(B,U)-77D(B), (11) 

where the term proportional to ly models the viscous ef- 
fect, the term proportional to rj models the Joule effect 
and F stands for the forcing. The linear operator D is 
defined as follows: 



T»{X) = {klxukl 



'2 2^2 7 



, Kjv Xj^) £ 



^Af 



(12) 



Now, the conservation laws must be enforced. Assum- 
ing incompressibihity, in the ideal limit and in absence of 
forcing (F, 1^,77 -^ 0), the model is expected to conserve 
the total energy £*°* — £^ +£^ , the cross helicity Ti"^ and 
the magnetic helicity H™. In terms of shell variables of 
the model, the energies and the cross heHcity are defined 
for the original MHD equation as 



1 






7T (U|U) , 



2 <^l^> 
(UIB). 



(13) 

(14) 
(15) 



The definition of the kinetic helicity and the mag- 
netic helicity requires the expressions for the vortic- 
ity O = (oi,...,OAr) and the magnetic potential vector 
A = (ai, ...,aAr). These quantities are not trivially de- 
fined in shell models since they require the use of the 
curl operator. Nevertheless, they should be linear func- 
tion of the velocity and magnetic field respectively. The 
kinetic helicity and the magnetic helicity are then defined 
as follows 



W" = (U|0) 

n^ = (AiB) 



(16) 

(17) 



In terms of conservation laws, the property[T]imposes the 
following constraints that correspond to the conservation 
of the kinetic energy and the kinetic helicity respectively 
by the first quadratic term in the U equation: 



(Q(U,U)|U) =OVU 
(Q(U,U)|0) =OVU 



(18) 
(19) 



Here, the notation "VU" must be understood as "for all 
possible values of the shell variables U as well as O that is 
defined by U". The conservation of the magnetic energy 
by the first quadratic term in the B equation (property[2]) 
imposes: 



The conservations of the total energy and of the cross 
heHcity respectively correspond to 

(Q(U,U)-Q(B,B)|U) + 

(W(U,B)-W(B,U)|B) =0 VU,B, (21) 
(Q(U,U)-Q(B,B)|B) + 

(W(U,B)- W(B,U)|U) =0 VU,B. (22) 

These two constraints are equivalent since the second 
is obtained simply from the first under the exchange 
(U,B) -^ (B,U). Hence, the general procedure adopted 
here shows that in the ideal Hmit, for a shell model with 
the structure l|10Hlip . the conservation of the total en- 
ergy f *°' implies the conservation of the cross helicity 
T-C^ and vice versa. Moreover, taking into account the 
constraints (fTSl and ([20l) . the conservation of the total 
energy and cross helicity reduces to: 

(Q(B,B)|U)-K(W(B,U)|B)=0 VU,B, (23) 

Finally, the conservation of the magnetic heHcity imposes 
the condition: 



(W(U,B)|A) + (W(B,U)|A) =0 VU,B, 



(24) 



Again, the notation "VU, B" must be understood as "for 
all possible values of the shell variables U and B as well 
as O and A that are defined by U and B respectively". 
The specific form of the nonlinear terms in the general 
shell model l|10mip can not be defined further without 
giving explicit definitions for O and A. The choice of the 
interactions retained in the nonlinear terms (for example: 
first neighbouring shell or distant shell interactions [la, 
NLTI I) must also be made expHcit in order to reach the final 
form of the shell model. An example will be treated in 
Section IVl 

If the shell model has to reproduce aH the symme- 
tries of the original MHD equation, the foHowing equality 
could also be imposed: 



W(X,X) ^Q(X,X) 



(25) 



(W(U,B)|B) =0 VU,B. 



(20) 



It is a consequence of the particular way of writing the 
MHD equations in which all nonlinear terms, including 
those appearing in the magnetic field equation, are made 
expHcitly divergence free through the appHcation of the 
projection operation ^. In the example treated in Sec- 
tion IIV^ this equality appears as a direct consequence 
of the other constraints imposed on the structure of the 
shell model. Nevertheless, if the present approach is ap- 
plied to more complex shell models for MHD, it might 
be interesting to keep the equality l|25p in mind in order 
to simplify the nonlinearities as much as possible. 



III. ENERGY FLUXES AND ENERGY 
EXCHANGES 

A. Evolution equations for the shell energies 

The kinetic and magnetic energies associated with the 
shell s„ are defined as e" = (U„|U„)/2 and e^ = 



(B„|B„)/2. The evolution equations for these quantities 
are easily obtained in the inviscid and unforced limit: 

dtel = K = (Q(U, U) - Q(B, B)|U„) , (26) 

dtel = t!: = (W(U, B) - W(B, U)|B„) . (27) 

The quantity T^ corresponds to the energy transferred 
into the velocity field in shell s„ and coming from either 
the velocity or the magnetic fields. Since the first term of 
equation l(26|) conserves the total kinetic energy (cf. I|18p ). 
it is identified as the rate of energy T^" fiowing from the 
complete velocity field into the velocity field in the n-th 
shell. The second term of fM)) must then account for the 
energy coming from the magnetic field (T^^), i.e., 



Tr = (Q(U,U)|U„) 
T,f = -(Q(B,B)|U„ 



(28) 
(29) 



Similarly, T^ corresponds to the energy transferred into 
the magnetic field in shell s„ and coming from either 
the velocity field or the magnetic field. The first term 
of equation l(27|) conserves the total magnetic energy (cf. 
(|20l) ). and is identified with the rate of energy fiowing 
from the complete magnetic field to the magnetic field of 
the n-th shell. The second term of Eq. (|27|) corresponds 
to the energy fiowing to the B„ shell from the complete 
velocity field, i.e.. 



rr = (w(u,B)|B„), 

r^ = -(W(B,U)|B„) 



(30) 
(31) 



With this notation, the evolution equations for e^ and 
ej^ become (with dissipative and forcing terms): 



dte: = rr + Tf - 2 V kl el + P/ ., 



dtet = T, 



bb 



T. 



bu 



2 77 /c^ e^ , 



(32) 
(33) 



where P/ — (F|U„) is the kinetic energy injection rate 
into the shell s„ due to the external forcing. 

It is also convenient to introduce the following decom- 
position of the vectors of shell variables: 

X,< = (a;i,2;2,...,a:._i,a;„0,...,0)eC^, (34) 

Xf = {0,Q, ...,G,Xi+i,Xi+2, ■■■,xn) G C^, (35) 

X = X<+X,>. (36) 

where i can take any value between 1 and N. The ki- 
netic energy contained in the vector U< is simply given 
by EV-" = (U<|U<)/2 = E;=ie,". The magnetic en- 
ergy contained in the vector B^ is defined similarly. The 
evolution of these quantities are easily derived from the 
relation l|32H33p : 

n n 

d,Er ^T.^r+T. '^f - ^- + p^"- ' (37) 






(38) 



where P/^ = (F|U^) is the injection rate of energy in 
U5 due to the forcing, and D± = v {Ti(\5)\\5^) and 



(D(U)|U<) and 



^nn = V 0^0^)\^n) cire the dissipative terms for t/,j 
and P,f respectively. 



B. Energy Fluxes 



< 



The nonlinear terms in the equations P7II38P corre- 
spond to the nonlinear energy fiuxes that enter or leave 
the sphere of radius /cqA". These fiuxes can be further 
specified. Indeed, the first sum in the right hand side 
of the equation (|37l) comes from the quadratic Q(U,U) 
term which conserves the total kinetic energy. Hence, 
this first sum must correspond to the kinetic energy fiux 
toU<: 



n^>(n) fromU, 



n 



U>( 



Et;"" = (q(u,u)|u< 



(39) 



j=i 



The anti-symmetry property for the fiuxes can be used 
to define the opposite transfer: Il^'^{n) = —II^^{n). It 
simply expresses that the energy gained by U^ due to the 
nonlinear interaction is equal and opposite to the energy 
lost by Uf . The magnetic energy fiuxes can be similarly 
defined as 



n 



B> 
B< 



N==E^f = (W(U,B) 



IB! 



(40) 



The "cross" fiuxes between the velocity and the magnetic 
field can also be defined systematically. The second sum 
in the right hand side of equation (|38l) corresponds to 
the fiux of energy from B< to U and readily leads to the 
following definitions: 



ng<H=E^"" = -(W(B,U)|B<), (41) 

N 

ng>H= E rf = (W(B,U)|B>). (42) 



j=n+l 



Since these terms are linear in U, each of them can easily 
be split into two contributions related to U^ and U^ 
respectively: 

nK>(n) 
nK>(n) 

The formula l|39H46p shows that the various fiuxes can be 
defined univocally, almost independently of the structure 
of the shell model as long as the terms conserving kinetic 
and magnetic energy have been identified. It must be 
stressed that, at this stage, the exact expressions for the 
nonlinear terms Q and W are not needed. 



= -n^<(n) = 


-(W(B,U<)|B<), 


(43) 


= -n^<(n) = 


-(W(B,U>)|B<), 


(44) 


= -n^>(n) = 


-(W(B,U<)|B>), 


(45) 


= -n^>(n) = 


-(W(B,U>)|B>). 


(46) 



C. Shell-to-shell energy exchanges 

The expression for some of the energy exchanges be- 
tween two shells may be derived from the above analysis. 
For instance, the quantity T^" has been identified as the 
energy flux from the entire velocity fleld to the magnetic 
field associated to the shell s„. The expansion ^ for U 
can be inserted into the term T!!^ and leads to: 



N 



r^-5]-(W(B,U, 



|B,„ 



N 

Erjibu 
nm 7 



(47) 



where each term in this sum can now be identified as the 
shell-to-shell energy exchange rate from the velocity field 
in the shell Sm to the magnetic field in the shell s„: 



T 



bu 



(W(B,U„0|B, 



(48) 



Similarly, by inserting the expansion ^ for B into the 
term T^ , it is possible to identify the shell-to-shell energy 
exchange rate from the magnetic field in the shell Sm to 
the magnetic field in the shell s„ as follows: 



Tr„ = (W(U,B„)|B, 



(49) 



Since the quantities T^^ are a shell-to-shell energy ex- 
change rate (the notation xy is referred to as general 
exchange and it can take values uu, ub, bu or bb), the 
following anti-symmetry property is to be satisfied: 



rpxy 



rpyx 



(50) 



It is worth mentioning that the present analysis does 
not lead to a simple definition of the shell-to-shell kinetic 
energy exchanges T^^. This is due to the presence of 
three velocity variables in the expression for the U-to-U 
transfers that prevents a simple identification of the ori- 
gin of the kinetic energy fiux. Nevertheless, considering 
the relation l(25|) . the quantity T^" l(28|) can be rewritten 
as follows: 



rr = (w(u,u)|u„). 



(51) 



and, by analogy with the expression (|49l) . it is reasonable 
to adopt the following definition: 



T™=(W(U,U„0|U„) 



The shell-to-shell energy exchanges give a more refined 
picture of the dynamics in the shell model than the fiuxes. 
It is thus expected that these fiuxes can be reconstructed 
from all the TjJ^^. The general formula are given by: 



N 



i^n+1 J— 1 



^$<<H-EE^:^^ 
1=1 j=i 

N N 

nEW= E E ^- 



(53) 

(54) 
(55) 



i— n+l J— n-t-1 



As a direct consequence of the property l(50|) . the same 
anti-symmetry property holds for the energy fiuxes. 

In the next section we will focus on a specific model 
adopted by Stepanov and Plunian |lj|. We will derive 
the formulas for the energy fiuxes and compute them 
numerically. 



IV. STUDY OF A GOY SHELL MODEL FOR 
MHD TURBULENCE 

A. Derivation of the shell model 

The results derived in Sections Ull and Hill are valid for 
any shell model for MHD that use only one complex num- 
ber per shell for each field (velocity and magnetic) and for 
which the properties (1-3) are satisfied. As long as the 
vorticity and the magnetic potential vector have not been 
defined explicitly, it is not possible to specify further the 
exact structure of the shell model, i.e. the structure of 
the nonlinear terms Q and W. In this section, we revisit 
the GOY-like shell model for MHD turbulence studied 
by Stepanov and Plunian [lj|, and apply the formalism 
discussed in Sections Ull and Uni to this model. The shell 
model is defined by the following expressions for the non- 
linear Q and W terms: 

qniX,X) = ifc„(ai<+ix*^2 + a2<_ia;*+i 

+ a3<_2<-i)> (56) 

Wn{X,Y) = ifc„(/?i<+iy;+2 + /32<-iy,*i+i 






PaV. 



n+l^n+2 



+ P5y*n-lXl+i + PQyl-2X*n-l) ■ (57) 

This shell model is fully determined if the following defi- 
nitions for the vorticity and the magnetic potential vector 
are also adopted: 



Oi = {-lyuiki , 
tti = {-lybi/ki 



(58) 
(59) 



(52) ft: 



Imposing the conditions derived in the previous section 
from the various conservation laws ifTSl [191 [20l [23l and 
lead to the following values of the parameters ai and 



A-1 



a2 



-Qfl- 



A2 



as 



1 



/3i = ai 



P5 



X^ + X + 1 
2A(A + 1) 
A^ - A - 1 
2A3(A + 1) 
A2 + A - 1 



-ai- 



P2 



X'-X 



ai 



2A2(A-|-1) 
A^ + A - 1 
2A(A-t-l) 

A^ + A + 1 



ai- 



2A2(A-1-1) "^^ '2A3(A + 1) 

As discussed at the end of Section [III this shell model 



also satisfies the constraint l|25l) . It is indeed easy to ver- 
ify that these parameters satisfy the following equalities: 

Pi+ P4=ai, f32+ P5= "2 and /S^ + Pe = as. 



In order to verify that the model derived here is exactly 
the same as the model discussed in [l2|, the dynamical 
system l|10mip can then be rewritten after a few algebraic 
manipulations as 

dtUn = ikn fpn(U, U) - _p„(B, B) j - vk^.Un + /„ (60) 
dt6„ = zA:„(w„(U,B)-t-„(B,U)) - i^klh^ (61) 

where 

Pn(X,X) — ai ( 2;„+ia;„-|-2 TY' -^n-l^n+l 

~ T3 ^n-2^n~lj J (62) 

Wn(X,Y) = [^11+11/71+2"^ ^n-lVn+l 

+ <-2yn-l) ■ (63) 



With these coefficients, the model (|60H6ip is clearly the 
same as the one derived by Plunian and Stepanov [ia |. 
Our interpretation of some of the shell-to-shell energy ex- 
changes and the energy fluxes derived in the Section IIIII 
and computed in the next section differs from those of 
Stepanov and Plunian [l2|- When we compare the two 
approaches carefully, we flnd that the velocity to velocity 
energy flux H^y^ and the total fluxes are the same for 
both the formalism, but other fluxes involving the mag- 
netic field are different. This is due to the fact that the 
complete function W is never computed in [I4I because 
the property [2] was not used explicitly in the derivation 
of the shell model. In particular, the part of the bilinear 
term that conserves the magnetic energy in the magnetic 
field equation was not identified. It was not needed to 
derive completely the model coefficient. However, this 
identification is needed if the energy fiuxes have to be 
defined unambiguously, which is the main objective of 
this work, but not of the approach developed by Plunian 
and Stepanov. 



In Fig. [T] we present some of the energy fiuxes as a 
function of wavenumber. We present the energy fiuxes 
n[/>("-), ni3>(^) and n^^(n) that are the rate of energy 
leaving a wavenumber sphere of radius fc„. In addition 
we also report the energy fiuxes n]j< 
defined as 



and n^;^ that are 



ttU< 



ttU< 
^^B< 



ttU< 
^^B> I 



n 



u< 

all 



n 



u< 
u> 



n 



u< 

B< 



n 



u< 

B> 



n 



u< 



n 



u< 

U> ' 



(64) 
(65) 



Finally, the total energy fiux n> leaving the sphere of 
wave vector smaller than fc„, from either the velocity or 
the magnetic field, is also presented on Figure [ij 
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Our energy fiuxes n> and 11^^ are in good agreement 
with the corresponding fiuxes reported by Stepanov and 
Plunian [I4I , which is expected since these fiuxes are com- 
puted by equivalent schemes in both these work. However 
the energy fiuxes from the velocity field to the magnetic 
field and vice versa do not match because these fiuxes are 
computed differently in these schemes. Our fiux formulas 
have the advantage of being defined by a more systematic 
method compared to that of Stepanov and Plunian |[l2| . 

Energy transfers 
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B. Numerical results 

In order to compute the energy fiuxes, we simulate the 
shell model ^l^ with v = IQ-^ and rj = 10"^ The 
magnetic Prandtl number is then Pm = vjr\ — 10^"^. 
The shells ratio is taken to be the golden mean: A = 
(1 -I- \/5)/2. We take the number of shells as iV = 36 and 
apply non-helical forcing to S4, S5 and se according to 
the scheme prescribed by Stepanov and Plunian [I4I . The 
energy injection rate e is 1. We evolve the shell model till 
our system reaches a steady-state, and then we compute 
the energy fiuxes by averaging over many time frames. 
Up to slight differences due to the time steps or the initial 
conditions, several steady state results obtained from our 
simulation very well reproduce the results presented by 
Plunian and Stepanov |ia] . 



Figure 1: Energy fluxes in function of the logarithm of fc„. 
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We also compute kinetic, magnetic, and total energy 
spectra in the steady-state. We observe that till fc ~ 10^ 
both the kinetic and magnetic energy show power law 
behaviour with -2/3 spectral exponent consistent with 
Kolmogorov's spectrum (Fig. [2]). After k ~ 10'', the 
magnetic energy decays exponentially due to the Joule 
dissipation, while the kinetic energy continues to exhibit 
power law behaviour with the same spectral exponent 
of -2/3. Stepanov and Plunian [I4I report two different 
spectral regimes at steady state. They report -1 spectral 
exponent at lower wavenumbers for both velocity and 
magnetic field. For higher wavenumbers, they report ex- 
ponential decay for the magnetic energy and -2/3 spectral 



In the following section we summarize our results. 



CONCLUSION 
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Figure 2: Kinetic and magnetic energy spectra z^ = 10 and 
Pm = 10-^ 



index for the kinetic energy. Our result differs from that 
of Stepanov and Plunian at lower waveumbers. We be- 
lieve this discrepancy is possibly due to a short range of 
wavenumbers that makes the determination of the expo- 
nent quite difficult. 

To understand the two spectral regimes better we fo- 
cus on the energy fluxes in these regimes. In the first 
regime (the lower wavenumber), the energy flux flgj]"^ 
and n> are approximately equal. The kinetic and mag- 
netic energy spectra show that both the velocity and the 
magnetic fleld are active in this regime. Hence it is rea- 
sonable to expect Kolmogorov-like energy spectrum for 
MHD turbulence according to the existing MHD tur- 
bulence phenomenologies (for review, see [lJ|). In the 
second regime (the larger wavenumbers), the magnetic 
fleld is damped heavily by Joule dissipation, while the 
velocity field evolves according to the nonlinear term 
Q(U,U) which corresponds to the convective nonlinear- 
ity u • V u. Hence the turbulence here is almost purely 
hydrodynamics except that C/^ shells still supply energy 
flux to the magnetic shells of the dissipative regime (small 
amplitudes). Therefore, it is reasonable to expect Kol- 
mogorov's energy spectrum for the velocity fleld in this 
regime as well but with a smaller energy flux. Hence, our 
energy spectra and energy fluxes shown in Figs. [Hand El 
are consistent since they exhibit almost the same slope 
but different amplitudes. 



A general derivation of shell models for MHD has 
been proposed. The conservation of the traditional ideal 
invariants of three-dimensional MHD turbulence is ex- 
pressed as general constraints that must be satisfled by 
the nonlinear terms in the shell model. The conservation 
of the kinetic helicity and kinetic energy by the hydrody- 
namic shell model in absence of magnetic fleld also leads 
to constraints on the nonlinearities. The similarity be- 
tween the original MHD equations and the shell model 
is pushed one step further by identifying one term in the 
magnetic fleld equation in the shell model that conserves 
the magnetic energy. It corresponds to the advection 
of magnetic fleld by the velocity in the MHD equations. 
This procedure is presented using a very general formal- 
ism which leads to a number of interesting results. 

We show that the conservation of the cross helicity 
and the conservation of the total energy are equivalent 
in shell models. This equivalence is a direct consequence 
of the symmetries of the MHD equations expressed by 
the general properties 1-3 presented in Section [TTl 

The expressions for the energy fluxes that are valid in- 
dependently of the speciflc structure of the nonhnear cou- 
plings between the shell variables have been derived. The 
knowledge of these fluxes is quite important when the 
shell models are used to explore dynamo regime. These 
expressions could even be used to derive shell models that 
would maximise or minimise certain energy transfers de- 
pending on the physics that has to be modelled. 

Also, expressions for the shell-to-shell energy ex- 
changes are derived. Like in the original MHD equa- 
tions, the energy exchange mechanisms in shell models 
unavoidably involve three degrees of freedom (triadic in- 
teraction) [i3, [3, M, H [la [U, [13, [ii- it is thus not 

obvious to derive expression for shell-to-shell energy ex- 
changes that are viewed as energy transfers between only 
two degrees of freedom. Nevertheless, the formalism pre- 
sented in Section lUl yields a very natural identiflcation of 
most of these energy exchanges. The only exception con- 
cerns the U-to-U energy exchanges. A simple expression 
is however also proposed for these quantities by analogy 
with the B-to-B energy exchanges. 

Another property of the formalism presented here is 
the clear separation between the treatment of the con- 
servation law and the assumptions that have to be made 
to deflne both the magnetic and for the kinetic helici- 
ties. Because these helicities involve quantities that are 
deflned using the curl operator, they are not very well 
adapted to shell models. It is thus quite appropriate to 
clearly present the expressions for the vorticity and the 
magnetic potential as additional assumptions required to 
fully specify the structure of the shell model. 

The procedure has been applied to a speciflc class of 
shell models based on flrst neighbour couplings, know as 



the GOY model. It has been shown that the general 
constraints naturally leads to the already derived GOY- 
MHD shell model |12|. However, the interpretation of 
the energy fluxes appears to be simpler in the present 
formalism. 

Several extensions to this work could be considered. 
Shell models using distant interactions between the shell 
variables [17| could be analysed using the same formal- 
ism. Also, despite the fact that the presentation has 
been made for shell models with one complex number per 
shell and per field (velocity and magnetic) , extending the 
present formalism to shell models with more degrees of 
freedom should be quite obvious. Finally, it would be 
interesting to explore other shell models based on alter- 
native definitions for both the vorticity and the magnetic 
potential. 
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